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Abstract: 

We have used the University of Tasmania Mt Pleasant 26m radio telescope to investigate the polarisa- 
tion characteristics of a sample of strong 6.7 GHz methanol masers, the first spectral line polarisation 
observations to be undertaken with this instrument. As part of this process we have developed a new 
technique for calibrating linear polarisation spectral line observations. This calibration method gives 
results consistent with more traditional techniques, but requires much less observing time on the tele- 
scope. We have made the first polarisation measurements of a number of 6.7 GHz methanol masers 
and find linear polarisation at levels of a few - 10% in most of the sources we observed, consistent with 
previous results. We also investigated the circular polarisation produced by Zeeman splitting in the 
6.7 GHz methanol maser G9. 62+0. 20 to get an estimate of the line of sight magnetic field strength of 
35±7 mG. 

Keywords: masers — polarisation 
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1 Introduction 



Interstellar masers from a variety of molecules, in- 
cluding OH, water and methanol have been detected 
in interstellar space. The geometry of the methanol 
molecule gives it a very rich rotational spectrum with 
more than a thousand transit ions at frequencie s less 
than 1 THz (see for example ICragg et al.l I2005T I. To 
date more than twenty of these transitions have been 
observed to exhibit maser emission towards high-mass 
star formation regions. Although many of these transi- 
tions are relatively weak and rare, the 6.7 GHz 5i-6 A+ 
transition is the second strongest interstellar maser 
transition (after the 22 GHz transition of water) and 
has been detect ed towards more than 900 regions withi n 
our Galaxy fsee lCaswell et al.ll2O10l : lGreen et alJl20ld . 
and references therein). This transition is exclusively 
detect ed towards young, h igh-mass star fo rmation re- 
gions l|Minier et all 120031 : IXu et all 120081 ). The role 
of magnetic fields in regulating star formation is still a 
matter of active debate, measurements of the polarisa- 
tion properties of interstellar masers have the potential 
to reveal details about the orientation and strength of 
the field at very high resolution for a few sel ect lines 
of si g ht within star formation regions (e.g. IPodsonl 
120081 ; IVlemmings et alj|20ict ISurcis et al.ll201ll ). Al- 
though a few objects have been studied in detail, to 
date there have been few studies which attempt to look 
at the polarisation properties of 6.7 GH z methanol 
masers in general (the exce ptions being IVlemmingsl 
120081 ; IVlemmings et al.ll201ll l. 

Here we present the results of a preliminary study 
of both the linear and circular polarisation proper- 
ties of 6.7 GHz methanol masers undertaken with the 
University of Tasmania 26m Mt Pleasant radio tele- 
scope. The maximum percentage linear polarisation 
detected in previous 6.7 GH z methanol maser stud- 
ies ra nges from a few to 10% (EUingscn 20021; IPodsonl 
2008), while the fraction of circular polarisation is even 



smaller (< 0.5%) (|Vlemmingsl I2008T ). Even though 
6.7 GHz methanol maser emission can be very strong 
(peak flux densities in excess of 1000 Jy), the majority 
of the emission in the majority of sources is signifi- 
cantly weaker. Accurate measurement of the polari- 
sation properties of the masers hence requires careful 
calibration of the instrument. The observations under- 
taken here represent the first spectral-line polarisation 
observations with the Mt Pleasant 26m telescope and 
in Section 2 below we outline in detail the calibra- 
tion process undertaken. Rather than solving all the 
polarisation properties of the masers simultaneously 
it was more efficient to used different methods to de- 
termine the linear and circular polarisation properties 
independently. The technique we applied to determine 
the linear polarisation (Stokes Q and U) assumes that 
there is no circularly polarized emission (Stokes V). 
Although this is not strictly true, the level of circular 
polarisation observed in 6.7 GHz methanol masers is 
sufficiently small that it has no effect on our ability to 
measure linear polarisation at present (other measure- 
ment uncertainties play a much greater role). 

There have been several other published methods 
for calibrating a single dish radio telescope for po- 
larisation observations. The approach used in this 
work for linear polari sation is similar to that used by 
ICenacchi et~aH l|2009h . The main difference is that a 
highly linearly polarised calibrator source is not needed, 
calibration can be performed with just unpolarised 
sources. This is achieved by only expanding to first 
order and by applying some simplifying assumptions 
made based on the properties of 6.7 GHz methanol 
masers. The use of unpolarised sourc es as calibrators 
is a general result of lHa makcr (2000), but this paper 
details why and how this is applicable to these obser- 
vations with the Mount Pleasant 26m. While some of 
the assumptions m ade in this work are known to be 
slightly inaccurate l|van Stratenll200l ). there are other 
sources of error that limit the accuracy of the results. 
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Other methods include observations of sources over a 
wide range of parallactic angles s u ch as the calibra tions 
performed bylHeiles et all <|200lh , Ijohnstonl <|2002h and 
Ivan Stratenl (120041 1. 

We have structured the paper in the following way. 
Section 2 outlines the theoretical basis of the method 
we have used to measure the linear polarisation prop- 
erties of the observed methanol masers. Differences in 
the telescope, receiver systems and spectrometer sys- 
tems at each telescope mean that there is no standard 
procedure for doing this and so we describe in detail 
the process used here. Section 3 summarises the ob- 
servations, and Section 4 describes the calibration ob- 
servations and processing required to apply the theory 
outline in Section 2. In Section 5 we discuss the re- 
sults of our linear polarisation observations towards a 
sample of strong 6.7 GHz class II methanol masers. In 
Section 6 we outline the procedures required to mea- 
sure circular polarisation in spectral line sources and 
in Section 7 we describe the results of our circular po- 
larisation observations. Section 8 presents the conclu- 
sions. 

2 Spectral line linear polari- 
sation measurement with the 
Mt Pleasant 26m 

In order to make linear polarisation measurements with 
any instrument it is necessary to relate the signal re- 
ceived by the processing hardware (in this case a digital 
autocorrelation spectrometer) with the signals emitted 
by the source of interest. While the general process 
through which this is done is common to all systems 
the specific details of the system (e.g. telescope mount, 
receiver characteristics) mean that in practice it is an 
instrument specific procedure. It is necessary to con- 
sider each section of the signal path from the emit- 
ting source through to the processing hardware and 
determine the effect of each section on the polarisa- 
tion properties of the signal. In this section we have 
undertaken that process for observations of 6.7 GHz 
methanol masers made with the University of Tasma- 
nia's 26m radio telescope at the Mt Pleasant obser- 
vatory. We have done this in the broader context of 
the mathematical framework of radio polarimetry for 
a single dish radio telescope, which (with appropri- 
ate modification) can be applied to observations with 
other receiver systems, or at other telescopes. The aim 
is to produce a 4x4 Mueller matrix in the Stokes frame 
that includes all the significant transformations that 
occur to the signal, including those due to the error in 
the telescope receiver feed. This will allow for calibra- 
tion of the telescope for polarisation observations. The 
starting p oint for this is the mathematical fram ework 
set out by lHamaker. Breeman. fc Sauld l|l996h which 
provides the connection between 2x2 Jones' matrices 
in the signal domain with 4x4 Mueller matrices in the 
coherency domain. This allows for the transforma- 
tions that occur to the observed signal as it propagates 
through space and the telescope system to be described 
by a series of relatively simple 2x2 Jones matrices with 



their corresponding 4x4 Mueller matrices. We consider 
each section of the signal path separately, starting from 
the emission of the radiation by the source through to 
the input signals to the processing hardware. Outlined 
in this section are these Mueller matrices for each sec- 
tion of the signal path and their combined product, 
the total Mueller matrix for the system. 

2.1 Faraday rotation 

As the emitted radiation from the source propagates 
towards the observer ionised material in the interstel- 
lar medium causes the linear polarisation to change 
position angle (Faraday rotation). The amount of ro- 
tation depends on the wavelength of the radiation and 
so for continuum sources this effect can be measured 
by observing the polarisation at a range of frequencies. 
However, this method cannot in general be used for 
spectral line emission as it occurs at a single narrow 
frequency range. In theory Faraday rotation can be 
measured for methanol masers, as they are known to 
produce cospatial emission from multip le transitions 
(|Menten et al.lll992l : iNorris et all 1 19981 ). If it is as- 
sumed that the polarized component of this emission 
is the same then observations of both transitions can 
be used to determine the rotation measure along the 
line of sight from the observed maser. Our original 
intention was to attempt this using observations of 
both the 6.7 and 12.2 GHz methanol transitions, how- 
ever, technical difficulties with the 12.2 GHz receiver at 
the time of the observations meant that we were only 
able to make polarisation observations of the 6.7 GHz 
methanol line and hence the Faraday rotation could 
not be determined. However, we do not expect this 
to significantly effect our results, as using models o f 
the electron distribution in the Galaxy iDodsonl (2008) 
found that the Faraday rotatation is +0.8° at 6.7 GHz 
for the source G339.88-l.26. While this will vary from 
source to source the amount of Faraday rotation at this 
frequency is small enough to be negligible compared to 
other sources of uncertainty in our polarisation mea- 
surements. 

2.2 Coparallactic angle rotation 

As the source being observed moves across the sky it 
also rotates relative to the telescope receiver. Most ra- 
dio telescopes have an azimuth-elevation mount, how- 
ever the Mt Pleasant telescope is an X-Y mount tele- 
scope (with the X-axis aligned north-south). For an 
X-Y mounted telescope this rotation is called coparal- 
lactic rotation (as opposed to parallatic rotation for an 
azimuth-elevation mount). From spherical trigonom- 
etry the coparallactic angle rotation can be shown to 
be 

6 P = 90° + arctan ( . (1) 
\sm(h) sm(fl) / 

Where 9 P is the coparallactic rotation angle, h is the 
hour angle of the source and 5 is the declination of the 
source. It should be noted that this is actually inde- 
pendent of the latitude of the observatory, because an 
X-Y telescope is aligned north-south and east-west. 
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Physical rotations, like the coparallactic angle, pro- 
duce a two times angle rotation in the Stokes coordi- 
nates Q and U. Hence coparallactic rotation for a na- 
tive circular receiver results in the following Mueller 
matrix 



1 




V o 





cos 28 p 
— sin 29 v 











cos26 p 



(2) 







1 



Note that the S superscript refers to the Stokes frame 
in the coherency domain. 

2.3 Receiver properties 

Firstly the receiver angle relative to the telescope axis 
has to be known. This is a physical property of the 
telescope and is constant unlike the time dependent 
coparallactic angle rotation. This rotation is described 
by the Mueller matrix Mj| od anglc which has the same 
form as Mf opara . 



The type of receiver feed also has to be accounted 
for. For a native circular receiver this is typically 
done in one of two ways; either by a transposition 
of the Stokes vector from e s = {/, Q, U, V} to e s = 
{/, V, Q, U} or by a change in the transformation from 
the instrumental frame (correlator outputs) to the Stokes 
frame. In order to use the canonical Stokes vector 
e s = {/, Q, U, V} the second approach will be adopted 
in this work. 




lHamaker fe Bregmanl {l996) give the transforma- 
tion from the Stokes frame to the right-left circular 
coherency frame to be 



(3) 



2.4 Imperfect receiver feed 

The error in the receiver has to be accounted for via 
the feed erro r Mueller matrix Mf rrnr . The frame work 
described by lHamaker~B rcgman . fc Saultl (|l996f ) ex- 
plains how to get from 2x2 Jones matrices representing 
transformations to the voltages to 4x4 Mueller matri- 
ces representing transformations to the Stokes param- 
eters. This allows for a large reduction in the num- 
ber of unknown parameters involved in the calibration 
process, depending on the Jones matrix used and the 
order of terms kept in the expansion. 



given by 



S _1 (J. 



1 J error ) S 



(5) 



Where ® is the outer matrix product, de f ined in Ap- 
pendix A of Hamaker. Br egman. fc Saultl l|l996h . 

Performing this expansion to first order in d pro- 
duces the following Mueller matrix in the Stokes frame 



/ 1 Xl + J/l X 2 + J/2 \ 

xi 1 j/i — Xi 

X2 + J/2 1 J/2 — X 2 

\ Xi — J/l X2 — J/2 1 



Where d x — Xi + 1x2 and d y = j/i + ij/2- 



/ 

(6) 



fourth row represent leakages to and from Stokes V. 
These terms are due to the difference between d x and 
d y and represent the amount of ellipti city present in 
the circular receiver (jHeiles et al.l f2001). These terms 
can be neglected in this case as Stokes V is expected 
to be insignificant due to the properties of the sources 
being observed. Without the out of phase couplings 
the error in the feed can be represented with just two 
parameters and the following Jones matrix 



1 d 
d 1 



(7) 



Where d = di+i(fe; di = |(ki + j/i) and <fe = \{x2 + J/2). 
This results in the following feed error matrix in the 
Stokes' domain 
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2di 


2d 2 





2di 


1 










2d 2 





1 
















1 



(8) 



2.5 System properties 

The telescope has an intrinsic gain and system phase 
for each of the channels. This can be represented by 
the following Jones matrix 



9i 




92 e 



(9) 



The general Jones feed error matrix used bv lHeiles et al.l 
is 

1 d x 



dy 



(4) 



Where d x and d y represent the complex leakage be- 
tween the two dipoles in the receiver due to misalign- 
ment or non-circularity. d y is the complex conjugate of 
d y . The corresponding Mueller matrix Mf rror is then 



Where gi and (72 are real and represent the voltage 
gains of the two polarisation channels. A noise diode 
connected to a short dipole, located inside the receiver, 
prior to the probes used to extract the incoming radia- 
tion is used to determine the system gains. The dipole 
attached to the noise diode produces linearly polarised 
radiation with roughly equal power received by both of 
the receiver probes. The strength of the noise diode is 
measured relative to sources of known flux density to 
determine the relative gains of each of the polarisation 
channels. 

The phase difference ip S ya in Equation [5] typically 
arises from differences in the electrical path length that 
the signals have to travel between the receiver and the 
correlator. This can be also dealt with by observing 
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the correlated signal injected by the noise diode. In the 
presence of this strong, correlated signal, a difference in 
electrical path length produces a linear variation with 
frequency of the phase of the complex cros s-correlation 
prod uct between the two polarisations (|Heiles et al.l 
2001). A linear trend line can be fitted to this across 
the bandpass and the result subtracted from the phase 
of the on source cross-correlation observations. For a 
native circular receiver this subtraction is a rotation 
between Stokes Q and U. It should be noted that this 
process introduces a small error due to the probe with 
the noise diode attached radiating into the feed, and 
hence the feed error could potentially change the effec- 
tive observed position angle of the probe. The emission 
from the noise diode probe is however highly linearly 
polarised, and so this error in position angle should be 
insignificant. 

The use of a correlated noise signal allows for the 
Mueller matrix M^, s arising from the Jones matrix 
Jsys to be accounted for before further analysis of the 
data is attempted. 



2.6 Overall signal path 



3 Observations 



Observations were conducted using the University of 
Tasmania's 26m radio telescope located at Mount Pleas- 
ant. This telescope has an X-Y mount, native circu- 
lar receivers at 6.7 GHz, a main beam FWHM of 7.5 
arcminutes at 6.7 GHz and a typical pointing accu- 
racy of 30 arcseconds. All of t he obser v ed so urces 
have accurate posit ions given by I Casw ell (2009) and 
ICaswell et ail (|2010h . These sources are compact and 
have negligible structure at single dish resolution, mean- 
ing that off-axis polarisation properties are insignifi- 
cant. 



A 2-bit digital auto-correlation spectrometer was 
used to correlate the data. For linear polarisation 
observations a 4 MHz bandwidth with 2048 spectral 
channels was used. This produced 4 outputs; the left 
and right hand circular auto-correlations as well as the 
real and imaginary parts of the left-right cross cor- 
relation. For circular polarisation observations a 4 
MHz bandwidth with 4096 spectral channels was used, 
this produced the left and right hand circular auto- 
correlations only. 



The overall signal path and its transformations on the 
incoming Stokes vector to produce the observed Stokes 
vector can be represented as follows 



*~obs ^"-sys " error " feed angle ■"■'■copara^'i! 



(10) 



The constant rotation due to the feed angle can be 
absorbed into the rotation for the coparallactic angle. 
The relative gains and and system phase can be ac- 
counted for to deal with M£, B . The type of feed has 
already been accounted for by specifying S. 



s 

^obs 



I s e s 

1 copara L 



(11) 



Here 



is the time dependent part of the total 



Mueller matrix of the system and Mf rror is the un- 
known time independent part representing the feed er- 
ror that needs to be determined via calibration. These 
can be combined into a single Mueller matrix M s de- 
scribing the system as a whole. 
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"obs 



M' 



M = 



V 



1 




K 2 





2di 


cos(2#p) 


sin(2fl p ) 





2d 2 


-sin(20 p ) 


cos(2#p) 














1 



(12) 
(13) 



Where K x = 2d x cos(20 p ) - 2d 2 sin(20 p ) and 
K 2 = 2di sin(2flp) + 2d 2 cos(26» p ). 

M represents the total transformations that occur 
to the incoming signal after taking into account the 
type of telescope feed and calibrating for the system 
phase. This description of the telescope system can 
now be used to determine what the unknown feed error 
is and calibrate the telescope for linear polarisation 
observations. 



The majority of the observations reported in this 
paper were carried out over the period 2010 Septem- 
ber 28-30. The observing strategy involved repeated 
observations of the strong 6.7 GHz methanol maser 
sources G323.74-0.26, G339.88-l.26 and G351.42+0.64 
over a range of hour angles. Observations of the con- 
tinuum source Virgo A for calibration purposes and of 
several other southern 6.7 GHz methanol masers were 
also undertaken. On 2011 February 12 further observa- 
tions were undertaken, repeating the calibration obser- 
vations of Virgo A and a few of the 6.7 GHz methanol 
maser sources observed in September 2010. In addi- 
tion we also made repeated calibration observations of 
the strong continuum source PKS B 1921-293 (another 
largely unpolarised source) to perform additional tests 
on the repeatability, accuracy and reliability of the cal- 
ibration procedures. 



Observations of the masers were made by first mak- 
ing a short off-source reference pointing with the noise 
diode on to allow correction for ip sys , followed by an- 
other longer off-source reference observation with the 
noise diode off, and then an on-source observation (with 
the noise diode off) of the same duration. The shape 
of the cumulative receiver and correlator bandpasses 
were corrected by subtracting the off-source reference 
observations from the on-source, rather than the more 
typical division method, because the cross-correlation 
product was not always greater than zero. 



A follow up circular polarisation observation of G9. 62+0. 20 
was made during August 2011, this consisted of a total 
of 4 hours on-source, 4 hours off-source; split into 10 
minute alternating scans. 
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4 Determining the feed error Stokes parameters as follows 



There is no single preferred method for measuring the 
feed error and hence obtaining the total Mueller matrix 
(Equation I13|) for th e observing system. One method 
which has been used ijHeiles et al.l200ll ; |Johnstorj2 002') 
involves observing a partially polarised source over a 
wide range of hour angles. This produces rotation 
in the observed source relative to the telescope feed 
due to parallactic (alt-az mount) or coparallactic (X- 
Y mount) rotation. Sine functions are then fitted to 
the fractional Stokes parameters of the source against 
the parallatic (or for an X-Y telescope, the coparallac- 
tic) angle. Fractional Stokes parameters are used due 
to there typically being variations in the total system 
gain with elevation. This method is especially useful 
when using native linear receivers and/or attempting 
to use more terms in the feed error matrix. The prob- 
lem with this method of calibration is that it requires 
observations over a large range of hour angles, requir- 
ing significant amounts of telescope time. 

Another method which can be used to determine 
the feed error matrix is to observe one or more un- 
polarised sources. This requires much less observing 
time, but since Stokes V is dominated by errors in rel- 
ative gain for a native circular receiver only two prod- 
ucts can be used, fractional Stokes Q and U. This pro- 
duces 2 equations and so cannot be applied to systems 
with more than 2 unknown parameters in the feed error 
matrix. It does however allow for independent calibra- 
tion of each spectral channel across the entire bandpass 
(in contrast to the previous method when applied to 
maser emission). 

In this section we outline how we used these two 
techniques to determine the total Mueller matrix (Equa- 
tion I13|l for our system. This is then used to deter- 
mine the feed error of the 6.7 GHz receiver on the Mt 
Pleasant 26m telescope. Observations of unpolarised 
sources were used as the primary calibration method 
in this work, however observations over a range of co- 
parallactic angles were also completed to compare the 
two methods. 



4.1 Partially polarised source 

When calibrating using a partially polarised maser source, 
the incoming circular polarisation and hence Stokes V 
is assumed to be zero, resulting in the following Stokes 
vector 



/ /in 



Uiu 





(14) 



After applying the total Mueller matrix for the system 
M s in Equation \W\ to this vector the observed frac- 
tional Stokes parameters can be expressed in terms 
of coparallactic angle, the feed error and the source 



9obs 



Qobs _ 2di/ in + cos(26lp)<2i n + siri(2<9j,)[/i n 



Ita + KiQm + K 2 U itt 



(15) 



^ohs — 



Uobs _ grfgjto - sin(26> p )Q in + cos(29 p )U in 



Ita + KiQta + KaUta 



(16) 



As the feed error matrix was only expanded to first 
order, the approximation can be made that J b s ~ Jin 
and hence the expressions reduce to 



q obs = 2di + Jqf n + u? n cos (20 p + $i) (17) 



u obB = 2d 2 + J qf n + uf n cos (26 p + $ 2 ) 



(18) 



Where $i and $2 are phase shifts depending on Qi n 
and U in . 

Observations of three strong (peak flux density > 
1000 Jy) 6.7 GHz methanol masers G323. 74-0.26, G339.f 
1.26 and G351.42+0.64 were completed over a range of 
coparallactic angles over the course of a day. For each 
source we averaged the emission over a velocity range 
covering the peak exhibiting the greatest linear polar- 
isation to produce a fractional Stokes q and u value. 
Non-linear least squares fitting was then used to fit 
a function of the form Y = yrcos(2# p + B) + C to the 
fractional Stokes q and u plots against coparallactic an- 
gle. While the source Stokes parameters are unknown 
at this stage, the A coefficient from this fit gives an 
estimate of the fractional linear polarisation present. 
The C coefficients from the fits are equal to the 2d\ 
and 2d,2 terms from the feed error matrix. The fits 
to the observations are shown in Figures [TJ{3] and the 
derived error matrix terms are given in Table [1] along 
with their formal uncertainties from the fits. The er- 
rors due to radiometer noise for each point in Figures 
[T}[3] are smaller than the size of the points used and 
hence are not shown. 



G323.74-0.26 fractional Stokes q and u 




-60 -40 -20 20 40 

Coparallactic angle (degrees) 



Figure 1: Fractional Stokes q and u against copar- 
allactic angle with fit for G323.74-0.26 



It should be noted that while the fit to fractional 
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G339.88-l.26 fractional Stokes q and u 




40 60 



Coparallactic angle (degrees) 



Figure 2: Fractional Stokes q and u against copar- 
allactic angle with fit for G339.88-l.26 



G351.42+0.64 fractional Stokes q and 




' -30 -20 -10 10 20 30 40 50 60 




-30 -20 -10 



10 20 30 40 50 

Coparallactic angle (degrees) 



Figure 3: Fractional Stokes q and u against copar- 
allactic angle with fit for G351.42+0.64 



Table 1: Coparallactic angle fit results 
2010 September 28-30 



Source 


2di 


2d 2 


G323. 74-0.26 
G339.88-l.26 
G351.42+0.64 


0.028±0.017 
0.030±0.012 
0.040±0.010 


0.156±0.006 
0.154±0.004 
0.126±0.002 


Mean 


0.033 ±0.007 


0.145 ±0.017 



The resulting observed fractional Stokes q and u are 
then simply given by 



<7obs — 



-/obs 



= 2di 



f^obs j 

^obs = -J = 

Jobs 



(20) 
(21) 



A continuum calibrator source, such as Virgo A, can 
then be observed to allow for calibration across the en- 
tire bandpass, from a single on-off telescope pointing. 
A linear trend line is fitted to the fractional Stokes q 
and u across the bandpass to reduce the effect of noise 
on the data, this trend line is then used to calibrate 
each channel individually. The means can work just 
as well however as the slope across the bandpass was 
typically very small. Observations of both Virgo A 
and PKS B 1921— 293 were made, with repeated ob- 
servations of PKS B 1921—293 undertaken to test the 
repeatability of the calibration method. Extragalactic 
radio sources such as Virgo A and PKS B 1921—293 
do sometimes exhibit linear polarisation at the level of 
a few percent in their core components. In single dish 
observations, the spatial blending within a single beam 
of these compact polarised components, is further di- 
luted when combined with the largely unpolarised kpc 
scale emission and together this reduces the fractional 
linear polarisation significantly. To first order these 
sources can be considered unpolarised to our system. 



Stokes q and u for G323. 74-0.26 in Figure [TJ does not 
appear as to be as good as the fits seen in Figures 
and [3] the errors in the fits are actually very similar. 
The error just appears large in Figure [1] due to lower 
fractional linear polarisation in G323.74-0.26 (around 
~ 2% compared to an average of ~ 6% for the other 
two sources). Another point to note is that while the 
range of coparallactic angles observed is less than 180° , 
the fits are restricted to a known angular frequency 
and so the length of the data set is sufficient to obtain 
accurate fits. 



4.2 Unpolarised source calibration 

When calibrating polarisation using an unpolarised source 
the incoming Stokes Q, U and V are assumed to all be 
zero. 

/ /in 

ein= I ] il'M 

V o 



Table 2: Unpolarised source calibration results 
2010 September 28-30 



Source 


2di 


2d 2 


Virgo 


0.0442±0.0005 


0.1383±0.0004 



2011 February 12 



Source 




2d 1 


2d 2 


Virgo A 




0.0795±0.0004 


0.0676±0.0004 


PKS B1921- 


293a 


0.093±0.001 


0.056±0.001 


PKS B1921- 


293b 


0.099±0.001 


0.065±0.001 


PKS B1921 


293c 


-0.177±0.001 


0.192±0.001 


Mean 




0.090 ±0.010 


0.063 ±0.006 



The errors quoted in Table [2] are the standard er- 
rors due to radiometer noise. The uncertainty in the 
mean is taken as the standard deviation of the values 
of 2di and 2^2 determined. The final result has an 
uncertainty significantly larger than that of the indi- 
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vidual observations, indicating that there is another 
unknown source of error besides radiometer noise af- 
fecting the results. The source of this error is possi- 
bly physical in nature, based on the direction the tele- 
scope is pointing. The Mount Pleasant 26m is an X-Y 
mount telescope, unlike alt-az mounts this results in 
differing orientations of the receiver platform relative 
to the ground and hence differing gravitational stresses 
on the receiver. These stresses could cause slight rota- 
tions or deformations of the receiver setup and slightly 
change the polarisation characteristics of the telescope. 

In Table [2] the a, b and c for the source PKS 
B 1921—293 refer to observations made at the begin- 
ning, middle and end of the observation run. The ob- 
servation at the end (labeled PKS B 1921— 293c in the 
table), produced very different values for d\ and d.2- It 
is not clear why this occurred and so this result has 
been included in the table, but was not included in 
the mean or further analysis. If it is due to the pre- 
viously mentioned physical effects then it would indi- 
cate they can be far larger than expected. Its presence 
does however indicate that multiple calibration obser- 
vations should be made so as to avoid possible individ- 
ual anomalous results. Another point of interest is the 
fact that the measured d\ and (fe changed significantly 
in between September 2010 and February 2011. Al- 
though the basic observing system and strategy were 
the same for the two sessions individual components 
in the system are likely to have been changed. Further 
work should be undertaken to monitor the magnitude 
and rate of change in the polarisation calibration of 
the system. 

As can be seen comparing the results for 2010 Septem- 
ber in Tables [1] and [2 there is fairly good agreement 
between the values determined for di and d,2 between 
the two methods. This indicates that the observation 
of unpolarised continuum sources is a valid calibration 
method for this type of telescope, though multiple cal- 
ibration observations should be made to ensure accu- 
rate and reliable results. 



Stokes parameters for G339.88-l.26 at 6.7 GHz 



5 Linear polarisation results and 
discussion 

The linear polarisation observations conducted are shown 
in the appendix. These are divided up into three groups; 
the sources which were observed in both September 
2010 and February 2011, the sources with significant 
linear polarisation observed only in September 2010 
and the sources observed in September 2010 which do 
not show significant linear polarisation. Significant lin- 
ear polarisation in this case is defined as greater than 
5 times the RMS noise level in the Stokes q and u 
observations, allowing for a meaningful plot of posi- 
tion angle and fractional linear polarisation. Tables [3] 
and [4] summarise these results for sources with signif- 
icant linear polarisation. Source coordinates in these 
tables a r e take n from the methanol maser catalogue in 
ICaswelll (|2009h . Example spectra for G339.88-l.26 are 
shown in Figures [3] and 




-38 -36 -34 -32 -30 -28 

LSRK Velocity (km/s) 



Figure 4: Stokes I, Q and U for G339.88-l.26, ob- 
servations made during September 2010 



Linear polarisation for G339.88-l.26 at 6.7 GHz 



-50 
-100 



-42 -40 -38 -36 -34 -32 -30 -28 




-38 -36 -34 -32 -30 -28 



Figure 5: Percentage linear polarisation and po- 
sition angle for G339.88-l.26, observations made 
during September 2010 



5.1 Comparison with other results 

Comparisons of several of the spectra were made with 
observations made using the ATCA in 1999 (Simon 
Ellingsen, Private Communication). Many of the fea- 
tures agreed to within one or two sigma for G9. 62+0. 20, 
G339.88-l.26, G345.01+1.79 and G351.42+0.64, which 
is a reasonable result considering the significant changes 
which have occurred in the Stokes I profile in the time 
interv al between the two sets of observations. iDodsonl 
(2008) also made polarisation observations of G339.88— 1.26 
at 6.7 GHz and these again are in approximate agree- 
ment with the observations made in this work. 

Observations were undertaken in February 2011 
to compare with the observations made in September 
2010. As can be seen from the appendix, these polari- 
sation spectra of the same source from the two sessions 
agree closely. This gives us confidence that the calibra- 
tion method is working correctly, even if the feed error 
properties of the telescope changed in between the two 
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sets of observations. 

5.2 Error estimation 

Three sources of error were included in our estimate 
of the total error in the Stokes Q and U values deter- 
mined. The first source of error considered is that due 
to the first order expansion in the feed error matrix. 
The first order expansion assumes that there is a small 
total leakage from Stokes I to Stokes Q and U. This 
was characterised by 

Ex = v /Q 2 + f/V(2rfi) 2 + (2d 2 ) 2 (22) 

The second source of error is that in the feed error 
matrix term. To estimate the uncertainty in each of 
these terms we used the distribution of values from the 
multiple independent measurements of these terms. 

E 2Q = JA(Mi) (23) 

E w = 7A(2d 2 ) (24) 

The final source of error we have considered is the mea- 
sured noise level in the Stokes Q and U spectra. This 
had to be included so that parts of the spectrum with 
no Stokes I signal did not confuse the final polarisation 
spectra. 

E 3 = RMS variation in Stokes Q and U (25) 

These errors were combined in quadrature to obtain 
the total estimated error in Stokes Q and U. 

AQ = ^E\ + E\ Q + El (26) 

AU = y/E* + E\ v + El (27) 

Uncertainty in the fractional linear polarisation and 
position angle were then determined by standard error 
propagation from the estimated Stokes Q and U errors. 

5.3 Discussion 

Recent wor k by lEUingsen ct al. (200?]); iBreen et al.l 
(|2010l . [2~01ll ) and others, has attempted to find an evo- 
lutionary timeline for the different maser sp ecies and 
transitions. In particular, IBreen et al ] (|2010h suggests 
that the more luminous 6.7 GHz methanol masers are 
associated with less dense dust clumps than the lower 
luminosity masers. The magnetic field is expected to 
scale proportional to the gas density and so we might 
expect to see some relationship between the maser lu- 
minosity and the polarisation properties (since these 
are thought to be produced by the magnetic field). Us- 
ing the kinematic distance to determine the luminosity 
of the 6.7 GHz masers we observed we looked to see if 
there is a relationship between the fractional linear po- 
larisation and the luminosity of the maser source. No 
such relationship was found for our sample of 6.7 GHz 
methanol masers. However, the primary selection cri- 
teria for the masers we targeted was that they have a 
peak flux density in excess of several hundred Jy (i.e. 
they are very bright masers). Further study involv- 
ing weaker, intrinsically less luminous masers should 



be undertaken to definitively test whether there is any 
relationship between the polarisation properties and 
the maser luminosity. 

The combination of comparisons with polarisation 
observations of 6.7 GHz methanol masers made with 
other telescopes and the checks for repeatability and 
self-consistency between different approaches of the Mt 
Pleasant observations, demonstrates that we are able 
to use this system to reliably measure the linear po- 
larisation of masers. The faster and easier calibra- 
tion method using unpolarised sources developed in 
this work appears consistent with the more commonly 
used approach of observing a partially polarised source 
over a wide range of hour angles. This development 
will hopefully allow further linear polarisation obser- 
vations of maser emission to be conducted with this 
telescope, or any other single dish radio telescope with 
native circular receivers. 

6 Circular polarisation 

Observations of circular polarisation with native circu- 
lar receivers are highly susceptible to errors in relative 
gain between the receivers. To deal with this the sys- 
tem gains have to be handled carefully. This is usually 
done by first assuming that any observed Stokes V is 
due to Zeeman splitting, typically due to a line-of-sight 
magnetic field. After making this assumption there are 
two main methods for determining the magnitude of 
this splitting covered in the literature, S-curve fitting 
and the running cross-correlation. 

6.1 S-curve fitting 

S-curve fitting involves fitting Stokes I and the deriva- 
tive of Stokes I (dl/dn) to the raw Stokes V spectra, 
which is the right hand circular autocorrelation minus 
the left hand circular autocorrelation. The form of the 
expression fitted is as follows 

V = al + b^- (28) 

Where a and b are the constants determined via the 
fitting. If the circular polarisation is due to Zeeman 
splitting then the part of Stokes V proportional to 
Stokes I is expected to be due to errors in the rela- 
tive gain, while the part of Stokes V proportional to 
the d erivative of Stokes I is due to the Zeema n split - 
ting. IVlemmings. Diamond fc van Langeveldel (2002) 
looked at synthetic Stokes V spectra and derived the 
following expression for magnetic field strength based 
on applying S-curve fitting to a Stokes V spectra. 

Vmax — Vmin _ 2A F _pi B\\ 
/max A«L 

Where Vmax and V m in are the maximum and mini- 
mum of the fitted Stokes V, Aul is the full width 
half-maximum (FWHM) of Stokes I and / max is the 
peak flux density of the maser feature. B\ | is the line 
of sight magnetic field strength and A F _ F > is the Zee- 
man splitting coefficient for the line observed. 
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This method produces a Stokes V spectrum as well 
as a fitted derivative of the Stokes I to the Stokes V, 
but requires the FWHM of the feature in Stokes I to 
determine the Zeeman splitting. This can prove prob- 
lematic in complex sources where the spectral peaks 
are not Gaussian and/or overlap in velocity; causing 
difficulty fitting them and blending resulting in reduc- 
tion in the magnitude of the Stokes V detected. 

6.2 Running cross-correlation 

Another method to determine the amount of Zeeman 
splitting directly without first forming a Stokes V pro- 
file is the running cross-correlation. The basics for per- 
forming a cross-correlation to determin e the amount 
of Zee man splitting are introduced by iModiaz et al.l 
(2005). The process involves forming three products 
between the left hand circular and right hand circu- 
lar autocorrelations at +1, and —1 velocity chan- 
nels of lag. As methanol is diamagnetic it exhibits 
weak Zeeman splitting and so the amount of splitting 
will typically be much less than the channel width. 
A quadratic can be fitted to these lags to determine 
how far in velocity the left and right hand autocorre- 
lations are shifted relative to each other. The Zeeman 
splitting coefficient of the line can then be used to de- 
termine the line of sight magnetic field strength. If x 
is the velocity shift between left and right hand circu- 
lar polar isations the line of si ght magnetic field is then 
given bv lModiaz et all l|2005h to be 



(30) 



V2A F _ F , 

IVlemming s ( 2008) extended the method of lModiaz et al 
(|2005h to use small velocity windows that are cross- 
correlated, these windows are moved across the spec- 
tra producing a running cross-correlation. This allows 
for variation in the magnetic field strength across the 
source to be determined, although smoothed depend - 
ing on the size of the windows used. IVlemmingsl ((20081 ) 
used 3 kms -1 wide windows, as noise considerations 
limited the use of smaller windows. In this work how- 
ever the only feature that provided a high enough sig- 
nal to noise ratio to perform this analysis on was the 
main peak of G9. 62+0. 20 after an 8 hour integration. 
A 1.6 kms" 1 wide window centered at the main peak 
was used, this avoided the emission from the secondary 
feature. 

6.3 Zeeman splitting coefficient 

The Zeeman splitting parameter for the 6.7 GHz methanol 
transition has not been determ ined exper imentally. Lab- 
oratory work has been done bv lJenl l|l95lft to determine 
the Lande-p factor of the 25 GHz methanol lines, this 
can be used to estimate th e splitt ing coefficient for the 
6.7 GHz line. IVlemmings! l|2008h gives an estimate of 
the splitting coefficient of the 6.7 GHz methanol line 
to be 0.049 km s _1 G _1 . It has recently e merged that 
an error was made in these calculations (Fish et alj 
1201 ll : IVlemmings et afll201lft and the correct result is 
an order of magnitude smaller. As a result of this we 
have repeated the calculations of the Zeeman splitting 



factor and o btain a value of 0.0048 km s ~ 1 G~ 1 (con- 
sistent with iFish et all ; IVlemmin gs et ahf). this is the 
value used in this work. Section 3 of IVlemmings et al.l 
(|201 lh gives a detailed discussion of the current state of 
knowledge about this important parameter. At present 
it would seem prudent to take the magnetic field esti- 
mates obtained in this and other work based on methanol 
maser observations as order of magnitude estimates, 
which may change if a better Zeeman splitting coeffi- 
cient is determined for the 6.7 GHz methanol line. 

6.4 Beam squint 

Offsets in the pointing between the left and right hand 
circular receivers can produce false Stokes V for ex- 
tended sources. Most emission from typical maser 
sources is however contai ned within a very small re- 
gion, less than 100 mas l|Caswelll Il997l : iMinier et al.l 
2000). This limits the false Stokes V due to beam 
squint to be typically much less than the true Stokes 
V, except when multiple sources are within the same 
beam. Comparisons with high resolution interferome- 
ter maps of the sources can help determine whether the 
Stokes V profile is due to beam squint. For example 
observations of G351.42+0.64 (NGC6334F) produced 
a profile consistent with significant beam squint due to 
strong emission in this source being present over an an- 
gular extent of several arcseconds. This is particularly 
an issue in this source as the emission from the two re- 
gions overlaps in veloc ity and is of comparable strength 
(jEllingsen et all 1 1996ft . G9.62+0.20 is also known to 
have two sites of emiss ion separated by around 13 arc- 
seconds fCaswcll 2009), however, in this case there is 
no overlap in velocity between the two sites and the 
omission at the offset site is much weaker. 

6.5 Other sources of error 

It should be noted that the approaches used to de- 
termine circular polarisation in this paper assume a 
perfect receiver. Non-circularity in the receiver feed 
will produce leakages from other Stokes products into 
Stokes V. Leakages from Stokes I are equivalent to er- 
rors in the relative gain and hence not an issue. Leak- 
ages from Stokes Q and U could potentially produce 
some false circular polarisation but this would typi- 
cally not be in a form mimicking the derivative of 
the Stokes I spectrum which might be mistaken for 
Zeeman splitting. The results of Section [5] also show 
that the magnitude of Stokes Q and U detected in 6.7 
GHz methanol masers would only produce measurable 
Stokes V if there were very high leakage. 

7 Circular polarisation results 
and discussion 

Both S-curve fitting and the cross-correlation method 
were applied to spectra from several sources. Most 
sources however did not produce a high enough signal 
to noise ratio in the observation time to give meaning- 
ful Zeeman splitting estimates. G9. 62+0. 20 was the 
only exception to this, with a line of sight magnetic 
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field strength of 35±7 mG via both the S-curve fitting 
and cross-correlation methods. The Stokes V spec- 
tra produced via S-curve fitting can be seen in Figure 
[6] Both of these methods were applied around the 
main peak only, avoiding the secondary emission. The 
Stokes V spectra for the secondary peak in Figure [6] is 
typical of beam-squint, which is expected as the emis- 
sion region is slightly off-set from t he main region by 
about 13 arcseconds l|Caswellll2069h . 



Circular polarisation in G9. 62 + 0. 20 at 6.7 GHz 




12 
L5RK Velocity (km/s) 



Figure 6: S-curve fitting for G9.62+0.20 



7.2 Comparison with other studies 

IVlemmings. Goedhart fc Gav lard (2009) looked at the 
Zeeman splitting in G9. 62+0. 20 and found an aver- 
age line of sight magnetic field strength of 11.0 ± 2.2 
mG for the main feature, with a similar result for the 
secondary feature. This was done however with the 
Zeeman splitting coefficient of 0.049 km s _1 G _1 . If 
the Zeeman splitting coefficient used in this work was 
applied then the field strength determined would be 
112±22 mG, much larger than the value obtained in 

this work. The spectral resolution used in this work 

was sl ightly higher than what was used by|Vlemmings, Goedh art fc 1 
(2009), by a factor of 1.25. This should have resulted 
in better sampling of the very narrow features of the 
S-curve spectra, resulting in an effective increase in de- 
tected magnetic field strength, rather than a decrease. 

The Zeeman splitting in G9. 62 +0. 20 is known to be 

variab le however, as discussed bv lVlemmings. Goedhart fc Gavlardl 
( 2009), so this may be the cause of the discrepancy be- 
tween the determined magnetic field strengths. 

7.3 Discussion 

While the circular polarisation results in this work 
are limited, they do at least show it is possible to 
make circular polarisation observations of methanol 
masers with small single dish telescopes. Further stud- 
ies of weaker masers using a 26m dish would involve 
many hours or days of observation for a single source. 
In practice large aperture single dish telescopes are 
required to make circular polarisation of all but the 
strongest methanol masers. 



7.1 Synthetic spectra simulations 

Monte Carlo simulations of synthetic spectra for G9.62+0 
were performed to estimate the error in the magnetic 
field strengths determined. These simulations involved 
first generating two Gaussians to represent the left and 
right hand circular autocorrelations. These had the 
same width and half the maximum of the Stokes I 
profile for G9. 62+0. 20 and were split in velocity by 
the shift determined via applying the cross-correlation 
method to the observed spectra. Simulated radiome- 
ter noise was then added to these, proportional to the 
system temperature in each channel based on the ob- 
servations. The sum of these was then used as the 
simulated Stokes I profile. The noise added was set so 
that the Stokes I baseline produced had the same noise 
level as observed in the Stokes I profile. The channel 
spacing in velocity was set to that of the observations, 
~0.044 kms -1 . The velocities of the channels relative 
to the generated spectra were also shifted randomly to 
reduce any effect under sampling may have had on the 
results. These synthetic spectra were then put through 
the S-curve fitting and cross-correlation methods to es- 
timate the magnetic field strength. This was repeated 
100 times and the standard deviation of the results was 
used as the error in the magnetic field estimates. 
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Table 3: Linear polarisation results, September 2010 



Source 


RA(2000) 


Dcc(2000) 


V pk (km s -1 ) 


Spfe (Jy) 


% Lin Polpfe 


Angle° fe 


Angle range 


G309.92+0.48 


13 50 41.8 


-61 35 10 


-59.7 


980 


5.9+1.3 


10+8 


-15-29 


G316. 64-0.09 


14 44 18.5 


-59 55 12 


-20.4 


120 


4.7+1.8 


51+8 


45-52 


G318. 95-0.20 


15 00 55.3 


-58 58 53 


-34.7 


560 


2.8+1.2 


20+14 


-33-22 


G322. 16+0.64 


15 18 34.6 


-56 38 25 


-63.0 


256 


2.5+1.6 


45+11 


43-60 


G323.74-0.26 


15 31 45.5 


-56 30 50 


-50.5 


2500 


2.3+1.5 


49+11 


-27-76 


G328. 81+0.63 


15 55 48.7 


-52 43 06 


-44.5 


280 


5.1+1.6 


-30+7 


-31-18 


G329.03-0.21 


16 00 31.8 


-53 12 50 


-37.0 


103 


4.1+1.8 


39+9 


34-39 


G335. 79+0.17 


16 29 47.3 


-48 15 52 


-47.7 


214 


5.6+1.3 


12+8 


5-33 


G339.88-l.26 


16 52 04.7 


-46 08 34 


-38.7 


1560 


0.6+1.2 


-23+56 


-32-67 


G345.01+1.79 


16 56 47.6 


-40 14 26 


-22.2 


306 


4.6+1.2 


14+10 


5-25 


G351.42+0.64 


17 20 53.4 


-35 47 01 


-10.4 


3430 


9.0+2.0 


48+5 


-89-73 


G353.41-0.36 


17 30 26.2 


-34 41 46 


-20.4 


84 


7.7+1.7 


-5+8 


-5-1 


G9.62+0.20 


18 06 14.7 


-20 31 32 


1.3 


5340 


2.3+1.2 


-23+15 


-23-11 




Table 4: Linear polarisation results, February 2011 






Source 


RA(2000) 


Dec(2000) 


V pk (km s -1 ) 


Spfe (Jy) 


% Lin Polp fc 


Angle^ 


Angle range 


G309.92+0.48 


13 50 41.8 


-61 35 10 


-59.7 


1090 


4.9+1.0 


19+6 


-16-44 


G323.74-0.26 


15 31 45.5 


-56 30 50 


-50.5 


3870 


3.4+0.8 


30+8 


-13-35 


G339.88-l.26 


16 52 04.7 


-46 08 34 


-38.7 


1660 


1.1+1.0 
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Polarisation Spectra 

Linear polarisation comparison spectra 

These sources were observed on both 2011 February 12 and 2010 September 28-30, to test the repeatability of the 
calibration method. The September observations for G339.88-l.26 appear in Figures [4] and [5] 

Stokes parameters for G339.88-l.26 at 6.7 GHz Linear polarisation for G339.88-l.26 at 6.7 GHz 





Stokes parameters for G309.92+0.48 at 6.7 GHz Linear polarisation for G309.92+0.48 at 6.7 GHz 
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Stokes parameters for G323. 74-0. 26 at 6.7 GHz Linear polarisation for G323.74-0.26 at 6.7 GHz 




Stokes parameters for G323.74-0.26 at 6.7 GHz Linear polarisation for G323.74-0.26 at 6.7 GHz 




Stokes parameters for G345.01+1.79 at 6.7 GHz Linear polarisation for G345.01+1.79 at 6.7 GHz 
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Stokes parameters for G351.42 + 0.64 at 6.7 GHz 



Linear polarisation for G351.42+0.64 at 6.7 GHz 
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Stokes parameters for G9. 62+0. 20 at 6.7 GHz 
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Linear polarisation for G9. 62 + 0. 20 at 6.7 GHz 
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Sources with significant linear polarisation 

These sources were observed 2010 September 28-30. Linear polarisation above the 5 times RMS level was considered 
to be a detection, plots of the position angle and fractional linear polarisation plots are restricted to those spectral 
channels which meet this criterion. 



Stokes parameters for G316.64-0.09 at 6.7 GHz Linear polarisation for G316. 64-0. 09 at 6.7 GHz 
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Stokes parameters for G328.81 + 0.63 at 6.7 GHz 




Linear polarisation for G328.81+0.63 at 6.7 GHz 
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Stokes parameters for G329.03-0.21 at 6.7 GHz 
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Stokes parameters for G335.79+0.17 at 6.7 GHz 
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Stokes parameters for G353.41-0.36 at 6.7 GHz 
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Weaker Sources 

The level of linear polarisation in these sources was not sufficiently high to allow position angle and fractional 
linear polarisation plots to be produced from the observations conducted 2010 September 28-30. 
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Stokes parameters for G336.43-0.26 at 6.7 GHz Stokes parameters for G337.61-0.06 at 6.7 GHz 
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